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Abstract 

In this paper we study the semiclassical limit for the pseudo-relativistic Hartree 
equation 

\J —£ 2 A + m 2 u + Vu = ( I a * \u\ p ) \u\ p ~ 2 u, in R N , 

where m > 0, 2 < p < V: M. N —► M is an external scalar potential, I a (x) = 

is a convolution kernel, CN, a is a positive constant and (N — l)p — N < 
a < N. For N = 3, a = p = 2, our equation becomes the pseudo-relativistic 
Hartree equation with Coulomb kernel. 
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1. Introduction 


In this paper we study the semiclassical limit (e —> 0 + ) 
relativistic Hartree equation 


dip 


= 


dt 


^\/—e 2 A + m 2 — nij ip + Vtp — j- * IPI 2 ^ V') 


for the pseudo- 
iel 3 (l) 


where ip: R x K 3 —» C is the wave field, m > 0 is a physical constant, e is the 
semiclassical parameter 0<eCl,a dimensionless scaled Planck constant (all 
other physical constant are rescaled to be 1), V is bounded external potential 
in R 3 . Here the pseudo-differential operator y/— e 2 A + to 2 is simply defined in 
Fourier variables by the symbol ^e 2 ^! 2 + to 2 (see [23]). 

Equation (1) has interesting applications in the quantum theory for large 
systems of self-interacting, relativistic bosons with mass m > 0. As recently 
shown by Elgart and Schlein [16], equation (1) emerges as the correct evolution 
equation for the mean-field dynamics of many-body quantum systems modelling 
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pseudo-relativistic boson stars in astrophysics. The external potential, V = 
V(x), accounts for gravitational fields from other stars. In what follows, we 
will assume that V is a smooth, bounded function (see [24, 19, 17, 18, 21, 
28]). The pseudo-relativistic Hartree equation can be also derived coupling 
together a pseudo-relativistic Schrodinger equation with a Poisson equation (see 
for instance [1, 32]), i.e. 

/ = (V — e 2 A + m 2 — to) ip + Vij) — Uil>, 


See also [14, 20] for recent developments for models involving pseudo-relativistic 
Bose gases. 

Solitary wave solutions ijj(t,x) = e ltx / £ u(x), A > 0 to equation (1) lead to 
solve the non local single equation 


V 


—£ 2 A + m 2 u + Vu 



( 2 ) 


where for simplicity we write V instead of V + (A — to). 

More generally, in this paper we will study the generalized pseudo-relativistic 
Hartree equation 


\/—e 2 A + m 2 u + Vu = ( I a * |uf) \u\ p 2 u, in M. N , (3) 

where to > 0, 2 < p < -AW, V : — > M is an external scalar potential, 


taix) = (x ^ 0), a G (0, N) 

is a convolution kernel and CN, a is a positive constant; for our purposes we can 
choose c/v,a = 1. For N = 3, a = p = 2, equation (3) becomes the pseudo- 
relativistic Hartree equation (2) with Coulomb kernel. 

We refer to [34, 9, 6, 30] for the semiclassical analysis of the non-relativistic 
Hartree equation. The study of the pseudo-relativistic Hartree equation (2) 
without external potential V starts in the pioneering paper [24] where Lieb and 
Yau, by minimization on the sphere {cj> £ L 2 (R 3 ) | f R3 \cj)\ 2 = M}, proved that 
a radially symmetric ground state exists in 17 1 / 2 (K 3 ) whenever M < M c , the 
so-called Chandrasekhar mass. Later Lenzmann proved in [22] that this ground 
state is unique (up to translations and phase change) provided that the mass M 
is sufficiently small; some results about the non-degeneracy of the ground state 
solution are also given. 

Successively, in [10] Coti-Zelati and Nolasco proved existence of a positive 
radially symmetric ground state solution for a pseudo-relativistic Hartree equa¬ 
tion without external potential V, involving a more general radially symmetric 
convolution kernel. See the recent paper [11] dealing existence of ground states 
with given fixed “mass-charge”. 

In [27] Melgaard and Zongo established that (2) has a sequence of radially 
symmetric solutions of higher and higher energy, assuming that V is radially 
symmetric potential. 
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The requirement that V has radial symmetry was dropped in the recent pa¬ 
per [8] , where a positive ground state solution for the pseudo-relativistic Hartree 
equation (3) is constructed under the assumption (N — 1 )p — N < a < N. 

To the best of our knowledge the study of the semiclassical limit for the 
pseudo-relativistic Hartree equation has been considered by Aki, Markowich 
and Sparber in [1]. Using Wigner trasformation techniques, they showed that 
its semiclassical limit yields the well known relativistic Vlasov-Poisson system. 

In the present paper we are interested to study the pseudo-relativistic Hartree 
equation in the semiclassical limit regime (0<eCl), using variational meth¬ 
ods. Replacing u(y) by £ 2 (i-p) u(ey), equation (3) becomes equivalent to follow¬ 
ing Hartree equation 

\J— A + m 2 u + V e (y)u = (/ Q * |u| p ) \u\ p ~ 2 u, in R w . (4) 

where V £ (y) = V{ey). In what follows we will assume that 

(V) V : R N -> R is a continuous and bounded function such that U m i n = 
inf R w V > —m and there exists a bounded open set O C R N with the 
property that 

Vo = inf V < min V. 

O do 


Let us define 

Jt = {y £ O | V(y) = V 0 } . 

We will establish the existence of a single-spike solution concentrating around 
a point close to . Precisely, our main result is the following. 

Theorem 1.1. Retain assumption (V) and assume that 2 < p < 2N/(N — 1) 
and (N—l)p—N < a < N. Then, for every sufficiently small e > 0, there exists 
a solution u e £ U 1 / 2 (] R N ) of equation (4) such that u E has a local maximum 
point y £ satisfying 

lim dist (sy £ ,^//) = 0, 

£->0 

and for which 

uffiy) < Ciexp(-C 2 \y - y e \) 

for suitable constants C\ > 0 and C 2 > 0. Moreover, for any sequence {e n } n 
with s n —^ 0, there exists a subsequence, still denoted by the same symbol, such 
that there exist a point yo £ with e n y £n —> yo, and a positive least-energy 
solution U £ lL 1 / 2 (R Ar ) of the equation 

V-A + m 2 U + V 0 U= (I a * U p ) UP- 1 

for which we have 

u s n ( V ) = U{y-y Sn )+ 3% n (y) ( 5 ) 

where lim„_> +oc \\& n \\ H i /2 = 0. 
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To prove the main result, we replace the nonlocal problem (3) in with a 
local Neumann problem in the half space R^ +1 as in [10] (see [4]). We will find 
critical points of the Euler functional associated to the local Neumann problem 
by means of a variational approach introduced in [2, 3] (see also [7]) for nonlinear 
Schrodinger equations and extended in [9] to deal with non-relativistic Hartree 
equations. 

In the present paper the presence of a pseudo-differential operator combined 
with a nonlocal term requires new ideas. As a first step, we need to perform a 
deep analysis of the local realization of the following limiting problem 

yf—K + m 2 u + au= ( I a * |u| p ) \u\ p ~ 2 u (6) 

with a > —m. This equation does not have a unique (up to translation) positive, 
ground state solution, apart from the case p = 2, N = 3. Nevertheless we can 
prove that the set of positive, ground state solutions to the local realization 
of equation (6) satisfies some compactness properties. This is the crucial tool 
for finding single-peak solutions which are close to a set of prescribed functions. 
Even if we use a purely variational approach, we will take into account the shape 
and the location of the expected solutions as in the reduction methods. 

Recently the existence of a spike-pattern solution for fractional nonlinear 
Schrodinger equation has been proved by Davila, del Pino and Wei in the 
semiclassical limit regime (see [15]). The authors perform a refined Lyapunov- 
Schmidt reduction, taking into advantage the fact that the limiting fractional 
problem has an unique, positive, radial, ground state solution, which is nonde¬ 
generate. 


Notation 

• We will use | • \ q for the norm in L q , and || • || for the norm in H 1 ( M^ +1 ). 

• Generic positive constants will be denoted by the (same) letter C. 

• The symbol R+ +1 denotes the half-space {(x,y) \ x > 0, y £ R N }. We 
will identify the boundary c>M+ +1 with R^. 

• The symbol * will denote the convolution of two functions. 

• For any subset A of R N and any g > 0, we set A e = {y | dist (y, A) < p}. 

• For any subset A of R N and any g > 0, we set A g = {y \ gy £ A}. 

2. Preliminaries and variational setting 

The realization of the operator \Jm 2 — e 2 A in Fourier variables seems not 
convenient for our purposes. Therefore, we prefer to make use of a local real¬ 
ization (see [10, 4]) by means of the Dirichlet-to-Neumann operator defined as 
follows. 
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For any e > 0, given u £ S(R N ), the Schwartz space of rapidly decaying 
smooth functions defined on R w , there exists one and only one function v £ 
<S(M^ +1 ) such that 


f —e 2 Av + m 2 v = 0 in R ^ +1 
\ v(Q,y) = u{y) for y £ R N = dR+ +1 . 


Setting 


TMv) 


we easily see that the problem 


dv 

£ dx 


( 0 , 2 /), 


f -£ 2 Aw + m 2 w = 0 in R ^ +1 

\ w(0, y ) = T e u(y) for y £ dR* +1 = {0} x R N ~ R N 


is solved by w(x, y) = y). From this we deduce that 

T e (T e u)(y ) = -£^(0,2/) =£ 2 |^(0,y) = (-£ 2 A y v + m 2 v) (0 ,y), 

and hence T e oT e = (—E 2 A y +m 2 ), namely T e is a square root of the Schrodinger 
operator —£ 2 A y + m 2 on = cM+ +1 . 


From the previous construction, we can replace the nonlocal problem (3) in 
R^ with the local Neumann problem in the half space R (( +1 

— e 2 Av(x, y) + m 2 v(x, y) = 0 in R ^ +1 

-£§y(°, 2 /) = ~V(y)v(0,y) + (I a * |u(0, -)\ p ) \v(0,y)\ p ~ 2 v(0,y) for y £ R N . 

Setting v e (x,y) = e 2 <~ 1 -p'> v(ex, sy) and V e {y) = V(sy), we are led to the local 
boundary-value problem 

( ~Av e + m 2 v e = 0 in K + +1 

l = ~V e (y)v e (0,y) + ( f a * K(0,-)l P ) K(0, y)\ p ~' V(0, y) for y £ R N 

We introduce the Sobolev space H = H 1 (R 1 ^ +1 ), and recall that there is a con¬ 
tinuous trace operator 7 : H —» iJ 1 / 2 ( R N ). Moreover, this operator is surjective 
and the inequality 

l7(«)lp < P\ v \ P 2{p-i) 

holds for every v £ H 1 (R^ +1 ): we refer to [33] for basic facts about the Sobolev 
space H 1 / 2 (R n ) and the properties of the trace operator. 

Reasoning as in [ 8 , Page 5] and taking the Hardy-Littlewood-Sobolev in¬ 
equality (see [25, Theorem 4.3]) into consideration, it follows easily that the 
functional S e : H —>• M defined by 

£ e {v) = 7 f |Vw| 2 dxdy + [ v 2 dx dy 

2 J r n +1 2 J R n+i 

+ \ ( V e {x)^{v) 2 dy - -t- f 

Z Jip J R N 


dv 

dx 
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(I a *h(v)n\l(v)\ p dy 





is of class C 1 , and its critical points are (weak) solutions to problem (4). 


3. Compactness properties for the limiting problem 

For a > —m, the equation 

\J — A + m 2 u + au= ( I a * |u| p ) \u\ p ~ 2 u (7) 

plays the role of a limiting problem for (4). Its Euler functional L a : H — > R is 
defined (via the local realization of Section 2) by 

L a (v) = 11 (|Vv| 2 + m 2 \v\ 2 ) dxdy 

+ ^ [ 17 {v)\ 2 dy-^-[ {I a * |T(^)r) \l{v)\ p dy. 

1 J R™ 1 V JW* 

We define the ground-state level 

m a = inf { L a (v) \ L' a (v ) = 0 ,v £ H \ {0}} 

and the set S a of elements v £ H \ {0} such that v > 0, L a {v) = m a , and for 
every x > 0 : 

max v(x, y) = v(x, 0 ). ( 8 ) 

ye 


Proposition 3.1. The set S a is non-empty for any a > — m. 

PROOF. The proof is indeed standard, and we will be sketchy. First of all, we 
invoke [11, Lemma 2.1] to deduce that ground states of L a correspond to ground 
states of the functional C a : H 1 / 2 (ML n ) —> M defined as 


£a{u) = i [ ( \J (m 2 - A ) 1/2 

z Jrn y 


— mu 


+ {a + m)\u\‘ 


[ (la* M p ) |«| p - (9) 

Z P J R« 


We claim that C a possesses a ground state. We fix a > —m and consider the 
minimization problem associated to (9) 

. f Irn | \/(m 2 - A ) 1 / 2 - tom | 2 + (a + m) |u | 2 

m a = mf —-j-. ( 10 ) 

uem,*(R»)\ { o} (f RN (I a *\ u \p)\u\p)r 


Since \fm 2 — A — m > 0 in the sense of functional calculus and a + m > 0, it 
follows easily that fh a > 0. As in [29, Proof of Proposition 2.2] we can show 
that rh a is attained. Since the quotient in (10) is homogeneous of degree zero, 
as in the local case we see that any minimizer of fh a is, up to a rescaling and 
a translation, a ground state for (9). Therefore the claim is proved, and in 
particular S a 7 ^ 0- It is easy to check that ground states are non-negative, and, 
as in [10, Theorem 5.1], actually strictly positive. 
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Remark 3.2. By [24, Formula (A.3)], the quotient to be minimized in (10) 
decreases under polarization. This implies, reasoning as in [29, Section 5] (see 
also [13]) that ground states are radially symmetric around a point of R N . 

For U £ S a , we write E a = L a (U). By an immediate extension of [31, Lemma 
3.17], the map a H > E a is strictly increasing and continuous. The following is 
the main result of this section. 

Proposition 3.3. The set S a is compact in H, and for some C > 0 and any a £ 
(—Vmin, m) D [ 0 , +oo) we have 

v(x, y) < Ce- (m - CT) V“ 2 + \y\ 2 e-™ (H) 

for every v £ S a ■ 

PROOF. If V £ S a , it follows easily from [10, Theorem 5.1] or [ 8 , Theorem 7.1] 
that v decays exponentially fast at infinity and (11) holds. Moreover, since 

m a = L a (v) = Q - (\Vv\l + m 2 \v\ 2 2 ) , 

S a is bounded in H. We claim that S a is also bounded in L°°(R+ +1 ). 

Indeed, by [10, Theorem 3.2] it follows that 'y(v) £ L 9 (R N ) for any q £ 
then also g(-) = —a'y(v) + ( I a * | 7 (u)| p ) \"l{v)\ p ~ 2 ^{y) £ L q (R N ) for q £ 
Following [5], we let u(x, y) = fg v(t, y) dt. It follows that u £ R 1 ((0, R ) 
for all R > 0. Arguing as in [10, Proposition 3.9], we can deduce that 
weak solution of the Diric.hlet problem 

( —Au + m 2 u = g in R ^ +1 
( u = 0 for y £ R^. 

where g(x,y) = g(y) for every x > 0 and y £ R w . We sketch the proof for 
the sake of completeness. Pick an arbitrary function rj £ C^°(R^ +1 ) and write 
u>t(x, y) = r/(x + t , y) for any t > 0. Then 


[2,oo], 
[2, ooj. 
x R n ) 
u is a 


( 12 ) 


/■+00 /■+ oo 

/ 0 Jo 


Vv(x, y) ■ V ??(x + t, y) dydx dt 

/» + 00 /* + oo p 


'0 J x 

/»+oo ps 


/R N 


Vv(x, y) ■ S7r](s, y) dy ds dx 


IQ Jo 


n \7v(x, y) ■ V? 7 (s, y) dy dx ds 

»JV 

v(x, y) dx\ ■ yrj{s, y) dy ds 


+oo 


0 J R N 


and this readily implies that 

/ (Vv-Vwt+iT^vwt) dxdy= / gwtdy. 

J r7 +1 J r n 
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An integration with respect to t from 0 to +oo gives 


X (Vu • V?? + m 2 ur) — grj) dxdy = 0, 

,N+1 

and hence the validity of (12) is proved. 

Moreover for any given R > 0 we can define u odd € H 1 ((—R,R) x R N ) and 
9 o M er)q>_ 2 L q ((-RR)*R N )hy 


^odd 


u(x,y) 

-u(-x,y) 


if x > 0 
if x < 0, 


g 0 dd{x,y) 


g{y ) if x > 0 

-g(y) if x < 0. 


It is easy to check as before that 

- Au odd + m 2 u odd = g odd in R N+1 . 


Since g 0 dd € L q ((—R,R) x R w ) for any q £ [2,+oo[, R > 0, we can invoke 
standard regularity results to conclude that 

Wodd G W 2 ’ q ((-R,R) xl f ) 

for every q > 2 and every R > 0, and hence u 0 dd € C' 1,/3 (R 7V+1 ), u £ C' 1, ^(R+ +1 ) 
and v = G C 0, ^(M^ +1 ) by Sobolev’s Embedding Theorem. Therefore g £ 
C°’P/(p~ i)(M 7V ) 1 and Schauder estimates yield u £ C 2 '^^ p ~ 1 ^ (K^ +1 ) and v £ 
C' 1 ,/ 3 /( p - 1 )(M^ +1 ). Moreover, the C 1 ’^- norm of v can be estimated by the L q - 
norm of g , which immediately implies that S a is a bounded subset of L°°(IR^ +1 ). 

Next, we claim that ^xn\(x,y)\^y+oo v (x,y) = 0 uniformly with respect to 
v £ S a - We assume by contradiction that this is false: there exist a number 
S > 0, a sequence of points (x n , y n ) £ M^ +1 and a sequence of elements v n £ S a 
such that x n + \y n \ +oo but v n (x n ,y n ) > 6 for every n. Let us write 
z n = (x n ,y n ), and call v n (z) = v n (z + z n ) for 2 = (x,y) £ R+ +1 . By the 
previous arguments, {v n }n is a bounded sequence in H nL 00 (M^ +1 ). Moreover, 
up to a subsequence, we can assume that v n —*■ v, v n v in H and locally 
uniformly in M^ +1 . As in [9, pag. 989], both v and v weakly solve (7). We now 
show that they are non-trivial weak solutions. The conclusion is obvious for v, 
since u n (0) = v n (z n ) > 5, so that 0(0) > 5. We consider instead v, and remark 
that [10, Eq. (3.16)] implies 


sup \v n (x,y)\ < C\i(v n )\2 e mx 

y£S. N 


for some universal constant C > 0. Hence S < v n (z n ) < | 7 (u n )|2 e~ mXrl , and 
the boundedness of j(v n ) in L 2 yields the boundedness of {x n } n in R. Without 
loss of generality, we can assume that x n —> x £ [0, +oo). Therefore, by (8), 

£ 

Vn(x, 0) > v n (x, y n ) > v n (x n , y n ) + o(l) > - 
by locally uniform convergence, and we conclude that v is also nontrivial. 


Now, for every n £ N, 


La{Vn) (2 2 p 

and 


/ (|Vw„| 2 +m 2 vl) dxdy + a / ^(v n ) 2 dy) = m a , 

Jm ) 


L a (v) > m Q , L a (y) > m a . 


If R > 0 satisfies 2 R < x n + \y n \, then 
m a = L a (v n ) 

> iiminf [ (| Vu n | 2 + m 2 u 2 ) 

\2 2pJ n-H-00 J B (q R ) 

+a / 7 {v n ) 2 dy 

J B(0,R)n({0}xR N ) 


2 1 ~- 2 - 2 ^ dxdy 


/ (IV^I 2 +m 2 v 2 ) 
\2 2p J ra-H-00 J B (o r.) 


+ ( — — 7-— ) liminf . 

^2 2p J n->+oo J B (0,R) 

+a / 7 {v n ) 2 dy 

J B(0,R)n({0}xR N ) 


2 , —2-2'! dxdy 


> 


1 1 


(| Vz>| 2 + m 2 v 2 ) dxdy + a / 7 (v) 2 dy 

J B(0,R)n({0}xR N ) j 


as 


.2 2py yjB(o,R ) 

= L a {v) + L a {v ) + o(l) = 2 m a + o(l) 

R —> +00. This contradiction proves that 

lim v{x,y) = 0 uniformly with respect to v £ S a . 
I (£,?/) 

From [10, page 70] it follows immediately that 

lim I a * |7(u)| p (y) = 0, uniformly w.r.t. v £ S a . 

|y|—H-00 

independent of v £ S a , such that \y\ > R a implies 
\I a *h(v)\ p (y)\h(v)(y)r 2 <^. 


(13) 


Pick R a > 0, ii 


As a consequence, 


f — Av + m 2 v = 0 in R+ +1 
I -&<-§« in {0} x {|j,| > Ra} 

As in [10, Theorem 5.1] or [8, Theorem 7.1], and recalling the uniform decay 
at infinity of (13), it follows that v decays exponentially fast at infinity, with 
constants that are uniform with respect to v £ S a . 
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We are ready to conclude: let {v n } n be a sequence from S a - Our previous 
arguments show that {v n } n converges — up to a subsequence — weakly to some 
v £ H, and this limit v is also a solution to equation (7). Fix 

r jv \ 

\ N(2 — p) + p J 

and split I a as I*+ 1%, where £ L r (R N ) and 7 2 £ L°°(R N ). This induces a 
decomposition of the non-local term <jY{v) = ^ / V 1 (v) + as 

,/V{v) = ^- f (I a *\'y(v)\ p )\'y(v)\ p dy 

2 P Jr jv 

JY 1 ^) = i [ (II * |7(V)D It (v)\ p dy 

2p J RN 

^ 2 {v) = ^- [ {II* ItMD \l{v)\ p dy. 

2p Jr w 

We obtain immediately that 


0 = lim 

n —^+oo 


/ (iVw^l 2 + m 2 v 2 ) dxdy - Jf(v n ) ] 

J r ” +1 J 

/ _ (|V,| 2 + m 2 v 2 ) dxdy — jV{v). 

J rT +1 


(14) 


We complete the proof by showing that lim ra _^ +00 jY{v n ) = ^Y(y). Now, by 
the Hardy-Littlewood-Sobolev inequality (see [25, Theorem 4.3]) 

|- - Jf 1 {V )| 

< f Il(x-y)\h{v n ){x)\ p \'y{v n ){y)\ p - \^(v){x)\ p \j(v){y)\ p \dxdy 


1 xl" 


1 Xl N 


Ia(x - y) |7K)(*)H7(«n)(v)l p - l7(«n)(*)| p |7(«)(l/)l p 


l7K)(®)H7(^)(y)l P - l7(v)(a;)H7(^)(2/)| P 


dx dy 


< 


1 xl" 


Ia( x - y)h(v n ){x)\ p ||t(w)(2/)| p - |tO)(?/)I p I dxdy 


+ 


/ R“xR» 


I l a (x - 2/)|7('y)(y)| P ||7WW| P - It(«)(®)I p I dxdy 


= 2 I a {x - y)\i{v n )(x)\ p ||7K)(y)| p - h(v)(y)\ p \dxdy 

J R N xI N 

< 2<7|£| r |7K)IV HtK)| p - ItMH - = 0(1), 

2r — 1 2r-l 

since |7(v n )| p — |t(v)| p —► 0 strongly in L^ -1 (R N ) by the choice of r. On the 
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other hand, 


\^ 2 {v n ) - JY 2 (v)\ 

< Halloo / ll7K)(a:)ri7K)(y)r - h(v){x)\ p \j(v)(y)\ p \ dxdy 
J 1" xl* 

and the conclusion follows as before. Since \im n ^ +00 ^V(v n ) = MT(v), equa¬ 
tion (14) yields limn^+oo ||u n || 2 = ||u|| 2 , and the proof is complete. 


4. The penalization scheme 


For 

S = ^ dist(^#, R n \ O) and /3 £ (0, 5) 

we fix a cut-off ip £ QJ°(M+ +1 ) such that 0 < p < 1 everywhere, <p(x,y) = 1 if 
x + \y\ < /?, and p(x,y) = 0 if x + |y| > 2/3. Setting p e (x,y) = < p{ex,ey), for 
any U £ >5V 0 and any point yo £ we define 


U'e° ( x , y) = p e (x, y - y ) U (x, y - y ) 
We also define, for all e > 0, 



if y e O e 
if y i o e 


and 

Qe(v) 


for v £ H. Finally, let 


Xel{v) 2 dy-1 


2p+l 

2 


+ 


r s {v) = S’eiv)+ Q e (v), V £ H. 

We want to find a solution, for e > 0 sufficiently small, near the set 

X e = { U*> | y 0 £ JJP and U £ S Vo } ■ 

We define the (trivial) path ip e (s) = sU for every s £ [0,1]. 

Lemma 4.1. There exists T > 0 such that T e (il> e (T)) < —2 for all £ sufficiently 
small. Moreover, 

lim max r e (-0 e (s)) = E v 

E-xO SG[0,T] 

where we recall that Ey 0 = Lv 0 {U) for U £ Syb¬ 


il 


PROOF. Indeed, by our definition of the penalization term Q E , by a simple 
change of variables and by the exponential decay of U at infinity, 


r e (^ e (s)) = S e {ip e {s)) 


|V *.(-)| 2 + ^ [ A(s ) 2 + ^[ V el {Us)f 

2 J R~ +1 2 

(4 * h ( A { s ))\ p ) |7(4(s))| p 

|VC/| 2 + n f f U 2 + \ f V 0l (U) 2 + 0(1)) s 2 
1 2 J rN + ! 2 J R N J 

[ (i a *\un\u\v+o(i)) s ^- 

Jwl n J 2 p 


2 J r" +1 
s 2p r 

2 P 

(l 

\2 4™+ 


where o(l) -> 0 as e -> 0 uniformly with respect to s. The conclusion follows 
easily. 


We are ready to introduce our mini-max scheme. For e > 0 sufficiently small, 
we define the set of paths 

<4 = W> e C([0,T\,H) I V’(O) = o, i>(T) = i!> e (T) = Tt/f} , 

where T > 0 is the number we found in Lemma 4.1. To this set we associate 
the min-max level 

C e = inf max r e 4(s)). 

V’S'i’e «e[o,r] 

By well-known arguments (see for instance [7, Proposition 3.2] for a proof in a 
local setting that extends smoothly to our case) it is possible to prove that 


lim C e = E v . 

v n u 


For a£l define the sublevel 


r“ = {w e h | r e («) < a}. 


Proposition 4.2. Let d > 0 be small enough, and let {£j}j be such that lim :) _ ) ._|_ 00 £j 
0 and let {u E; .} C Xf. be such that 

lim T e (v E )<E Vo , lim T' (v e ) = 0. 

j —>-+oo j—>-+oo 3 

Then there exist — up to a subsequence - {yj}j C R^, a point y G and 
U G Sy 0 such that 

lim |ejyj -y |=0 

J —>-+oo 

lim ||u £j - Vj)U{- ,• - yj) || = 0. 
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PROOF. In the proof we will drop the index j and write e instead of £j for 
simplicity. By Proposition 3.3, there exist Z £ Sy 0 . {y e } C .y# 3 and y £ 
such that y e —>• y as e —> 0 and 


( 

Ve\ ry ( 

Ve\ 

Ve - Pe y- 


■-J) 


< 2d for every e<1. (15) 


We set 


Vl,e = Pe (■, ■ - y ) Z (•, • - y ) , V 2 ,e =v e - V^ e . 


We claim that 

r E (w £ ) > r e (ui )£ ) + r e (^2, E ) + o(e). 

Suppose that there exist R > 0 and points 

'ife 2/3\ \ n (P_' 


y e &B\°-r—)\B 

£ £ J \ £ £ 


such that 


liminf / 7 (u e ) 2 dy > 0. 

E ^° JB(y e ,R) 

Set v e (x, y) = v e (x, y + y £ ) so that 


liminf / 7(u £ ) 2 dy > 0. 

e_v0 Jb(o.r) 


Ib(o,r) 

Up to subsequences, we can assume that 


(16) 


(17) 


lim eye = yo G B(y,2P)\B(y,P). 
£—>■0 


The sequence {u E } is bounded in id and hence in every L q (M. N ) with q < 
2N/{N — 1). As a consequence, v e — >• W weakly in H and strongly in L 1 9 oc (M JV ) 
for every q < 2N/(N — 1). By (17), W 0. Moreover, 

y/-A + m 2 W + V(y 0 )W = (I a * |W| P ) |W| P " 2 W. 

Choosing i? 1, 

liminf [ (|Vt; £ | 2 + m 2 v^) dxdy> ^ / (|VW| 2 + m 2 W 2 ) dxdy. 

e_>0 J(0,+oo)xB(y e ,R) 2 Jr«+! 

Since E a > Ei, whenever a > b, we have 


L v( yo )(W) > E V ( y0 ) > £V 0 . 
Hence, for some absolute constant cq > 0, 


liminf / (IV« £ | 2 + m 2 v 2 e ) dxdy > c 0 • L v{yo) (W) > c 0 • E Vo > 0, 

e ^° J(0,+oc)xB(y e ,R) 
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and this is a contradiction to the exponential decay at infinity of Z and the fact 
that 2/0 7 k V- 

Since such a sequence {y e } cannot exist, a Lemma of P.-L. Lions (see [26, 
Lemma 1.1]) implies that 

limsup / |7(v e )| N ~ 1 dy = 0. 

e-H) J 

This, the boundedness of {7(u E )} in L 2 and the Hardy-Littlewood-Sobolev in¬ 
equality imply 

lhl n ( f ( Ja * h( V e)\ P )h( v e)\ P dy 

- [ (I a *\j{vi,e)\ p )h(vi,e)\ p dy - [ (I a * \'y(v 2 ,s)\ p )\l(v 2 , e )\ p dy] = 0. 
J R« Jr n J 

If we write 


r £ (u E ) = r e (wi, e ) + r e (u 2 , e ) 


+ <Pe [X,y- —) (l - We (x,y - 

J(0,+oo)x(B(^,^)\B(^-,|)) V e/V V 

+ [ V e 'y(ip e (x,y-—))(l-'y(<p e (x,y- 

J B (Me. (Me. ? £) V V £ / / \ V V 

( T a * \l(Vs)\ P )\l(Ve)\ P dy 

2 P Jrn 

+ ^~f ( I a*\'y(vi,e)\ p )\'y(vi >e )\ p dy 

2 P J R« 

+ ^~ [ (da * \'y(v2,e)\ P )\l(v2, s )\ P dy + o(l) 

2 p J RJV 


y)) |Vu e | 2 died// 
y))) l(v s ) 2 dy 


as e —> 0, we deduce that (16) holds true. We now estimate r e (u 2iE ). There 
results 


r e (u 2jE ) > £ e (v 2 ,c) 

= \ [ |Vu 2jE | 2 dxdy + ^ f 147(^2 ,e) 2 dy 
2 Jr** 1 2 J r n 

- [ (da * \l(V2,e)\ P )\l(v2,e)\ P dy. ( 18 ) 

2 p J RJV 

For some constant C > 0 and using again the boundedness of {7(1^, E )} in L 2 , 

[ (da * \l(v2,e)\ P )\l(l%e)\ P dy < C'||u 2iE [|, 

J R N 

Now ( 15 ) implies that ||i; 2iE || < 4d for small values of s. Taking d > 0 sufficiently 
small uniformly with respect to e, we have 

JlKsf - f (da* h(v2,e)\ P )\l(v2, E )\ P dy > J||u 2 , E || 2 . 

Z Z P J R» O 
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Since the functional § e is uniformly bounded in X'J- for small e > 0, the pe¬ 
nalization term Q e is uniformly bounded in for small e > 0 as well. As a 
consequence, for an absolute constant C > 0, 


/i«\o £ 


7(^2 ,ef dy < Ce», 


and (18-19) imply r(u2, E ) > o(l) as e —> 0. 
Let us introduce 


vl, e (x,y) 


vi, E {x,y) 

0 


if y e O e 
otherwise. 


(19) 


For 2 U e (x,y) = v\ e {x,y + y E /e ), we can proceed as before and conclude that 
2D E converges weakly in L 9 (R^ +1 ), q < 2N/(N — 1), to a solution 2U of 

yZ-A + m 2 W + V(y)W = (/« * |QIJ| p )|Qir| p - 2 2n. 


We claim that 23J E converges to 2U strongly in H. As before, assume the 
existence of a radius R > 0 and of a sequence {z E } C such that z e £ 
B(y E /e,2p/e), 


liminf|z E — e 1 y e 1=0 and liminf / \l( v \ E )| 2 dy > 0. 

e_> ' 0 E_>0 Jb{z c ,r) 

Without loss of generality, ez E — > z £ O as e —> 0. Then 2B E (:r, y) = 2IJ E (:r, y + 
z E ) converges weakly in L q (R^ +1 ), q < 2N/(N — 1), to some 2U £ H that solves 

V 7 -A + to 2 2IJ + V(z)W = (I a * \W\ P )\W\ P ~ 2 W. 


and we obtain a contradiction as before. Again, 

lim if (I a *\ 1 m)\P)hm)\Pdy= f (I a * |7(2U)| P )|7(2IT)| P dy. (20) 

E_>U il» J RJV 

Hence 


limsupr e (v^ E ) > liminf r e (uj e ) 
e-vo ’ s ~>° 

> liminf ^ [ |V2B E | 2 dxdy 

£->0 2 J(o,+oo)xB(0,R) 

+ \ [ V(ey + y E )|7(2H e )| 2 dy 
A Jb(o,r) 

-f/ (i a *MW E )\ p )hm)\ p dy 

2 p J RJV 

>\[ |V2B E | 2 dxdy 

2 J(0,+oo)xB(0,R) 

+ \ v {y) [ hi^e^dy 

A Jb(o,r) 

-f/ (I a *hm\ p )hm\ p dy. 

2 P J R« 
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Since R > 0 is arbitrary, 


limsupr e (rij e ) > 
£—>-0 



|V2U| 2 dxdy+^-V{y)[ |'y(22J )| 2 dy 
z J e n 


(i a *\ 7 (w)nhm\ p dy 


= L V(9) m 

> E Vo . (21) 


Recalling (16), we find 

limsup (r E (w 2 , e ) +r E (r>J £ )) = limsup (r e (z; 2 , £ ) +r E (v iiE )) < limsupr E (i; E ) < E Vo . 

e— y 0 £—>-0 £—>-0 

Now r £ (M 2jE ) > o(l) yields 

lim n r E («l ) = e Vo . 

£—>•0 ’ 

What we have just proved entails that L v ^{ 2B) = Ey 01 and then y £ . As 

a consequence, 2U is, up to a translation in the y- variable, an element of Sy 0 , 
namely W(x , y) = U(x,y — 3 ) for some U £ Sy 0 and some 3 £ R N . 

Recalling that V >V(y) on the subset O and using the identity L v (y) (2R) = 

E Vo we get 

I \V<m\ 2 dxdy + V 0 f \ 1 m\ 2 dy-2p [ (I a * | 7 (2R)n| 7 (2IJ)| p dy 

J r ” +1 Jr n Jr n 

> limsup f |V2R £ | 2 dxdy + f V(ey + y e )\^){W e )\ 2 dy 

£->o Jr ™ +1 Jr n 

-2 pf (/ ct *|7(2ne)| p )l7(2R,)l P rfy 

Jr n 

> limsup f \\7W e \ 2 dxdy + [ V(y)\~/(W e )\ 2 dy 

e->o Jr ™ +1 Jrn 

-2 p[ (/ Q *|7(2n e )l P )l7(2Re)l p dy 

Jr n 

>[ \VW\ 2 dxdy + V 0 [ |7(2R)| 2 d2/-2p / (J Q *|7(2R)ni7 (2R)IW) 

J r " +1 Jrn Jrn 

and therefore 

lim / V(ey + y e )\"f(W s )\ 2 dy = [ V(y)\j(W)\ 2 dy. 

Jrn Jrn 

Using again the fact that V > V(y) on the subset O we conclude that 7 ( 2 R E ) —> 
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7(211) strongly in L 2 (M. N ). Finally, from (20), (21) and (22) we see that 


[ |V2IT| 2 dxdy+ [ V(y)\im\ 2 dy 

Jr ™ +1 Js. n 

> limsup f |V2h E | 2 dxdy + I V(y)\j(W e )\ 2 dy. 
e->0 JR ^ +1 J RJV 

The strong convergence of 2U E to 2U in H is now proved. Thus 

v{, E (x,y) = u(x,y-^~ 3) +o(l), 

and straightforward algebraic manipulations show that 

vi, e (x, v) = <Pe(x,y-‘!j-l)u(x,y-‘^-i)+ o(l) 

strongly in H. But Ey 0 > lim e _ >0 T E (u E ) and lim E _>o T e (ui iE ) = Ey 0 , so that 
lim £ _>.o T e (u 2 ,e) = 0 by (16). Using (18) and (19) we discover that V 2 , £ —t 0 
strongly in H. This completes the proof. 

5. Critical points of the penalized functional 

We are now ready to show that the penalized functional T E possesses a 
critical point for every e > 0 sufficiently small. 

Lemma 5.1. For d > 0 sufficiently small, there exist positive constants e 0 and 
u such that |T^.(u)| > w for every v £ T® e n (xf \ X^ 2 ^j and e £ (0, £0). 

PROOF. If not, for d > 0 so small that Proposition 4.2 applies, there exist 
sequences {£j}j with linij Ej = 0 and {v Ej }j with v e . £ X \ X £ . 2 satisfying 

lim r (w ) < E Vo and lim r' (v e ) = 0. 

J-t-j-OO J —>-+00 3 

Hence Proposition 4.2 applies and provides points y £j £ R N , y £ -jM and a 
ground state U £ Sy 0 such that 

lim |Ejyj -t /|=0 

J —>-+oo 

iim \\v E - <p E (•,• - yj)U(-,- - yj)\\ = 0. (23) 

j—>-+oo 

The definition of X e . implies lim^+oo dist (v e , X e .) = 0, and this contradicts 
the assumption v e . X E ( 2 . 

Let now d > 0 be chosen so that Lemma 5.1 applies. 

Proposition 5.2. For e > 0 sufficiently small, the functional P E has a critical 
point v e £ X E n Tf. 
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PROOF. Pick R 0 > 0 so large that O C ({0} x R N ) n B(0,R 0 ) and ^> E (s) £ 
Hq (B(0, R/e)) for any s £ [0, T], R > R 0 and e > 0 sufficiently small. We write 
D e = maxo< s <Tr £ (!/i £ (s)). By Lemma 4.1, there exists a £ (0 ,Ev 0 ) such that, 
for sufficiently small e > 0, 

T s (t/j e (s)) > D E - a implies ip e (s) £ X/ /2 n Hq (5(0, R/e)). 

We claim that, for sufficiently small e > 0 and R > Rq, there is a sequence 
{v^}„ C X d J 2 n rf* n H\{B{ 0, 4?/e)) such that r' e (v%) -)• 0 is ^(B(0, i?/e)) 
as n —> +oo. 

Arguing by contradiction, we assume that for sufficiently small e > 0 there 
exists a number aj{(e) > 0 such that 

|r e (t;)|>Ofl( e ) 

on X d/2 rr°* HHq (B(0, R/e)). With a slight abuse of notation, we will identify 
any v £ Hq(B(0,R/e)) with its extension to H as the null function outside 
B( 0, R/e). Applying Lemma 5.1, we find a number u > 0, independent of e > 0, 
such that |r' £ (u)| > w for v £ T// e D (Xjf \ X^ 2 ). By a classical deformation 
argument that starts from ib e , there exist some u £ (0, a) and a path ib £ 
C([0,T\,H) satisfying 

V’(s) = ip e (s) for tfj e (s) £ rf e_ “, ip{s) £ X d for ip s (s) £ rf e_tl 

and 

T s (ip(s)) <D e - y for every s £ [0,T], (24) 

Let C £ Co°(K^ +1 ) be a cut-off function such that ((x,y) = 1 for 0 < x < 8 
and y £ O s , C,(x,y) = 0 for x > 28 and y ^ O 26 , £(•,•) € [0,1], and |V£| < 2/5. 
For i/j(s) £ X// we denote tpi(s) = CeV’( s ) and ^(s) = (1 — ( e )i/j(s), where 
( e (x, V ) = C( £X i £ U )• We remark that we understand the dependency on e in the 
notation of if)\ and 1 / 2 - Observe that 

TeWs)) = r e (V>l (s)) + r e (i> 2 (s)) + Qe(l/(s)) ~ Q e (lf>l(s)) ~ Q E (lp 2 (s)) 

-^i (i a *hWs))\ p )hWsW 

2 V J R« 

+ -rz[ ( 7 « * \i{Ms))\ p )hi^i(s))\ p 
Z P J R N 

+ ^[ (i a *h(Ms))\ p MMs))\ p . 

2 V Jvl n 

The elementary inequality (h + k — 1)+ > (h — 1)+ + (k — 1)+ valid for h > 0 
and k > 0 immediately implies that 

Qe{i/>(s)) > Q e (V>i(s)) + Q e (ip 2 (s)) 

and, similarly to (19), we find that 

[ \l(^{s))\ 2 dy < Ce 6/p . (25) 

JR N \O e 
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On the other hand, writing k = (/ a *|7(V’(s))| p )|7(V’(s))| p -(/ a *|7(^i(s))| p )|7(V’i(s))| p - 
(/ Q * |7(V’2(s))| p )|7(i/’2(s))| p for simplicity, 

[ K = 2 f (I a * |7(V’(s))| p )|7(V’(s))| p 

Jr n J0™x(R n \0™) 

-2 [ {I a * |7(^ , (s))| p )|7(V , (s))| p 

7( 0 M\ 0 |)x(RJV\ 0 J) 

and from (25) via interpolation we deduce that 

lim f (I a * |7(V’(s))I p )|7(V’(s))| p = 0 (26) 

e ->°J 0* s x(R N \0* 8 ) 

lim/ (I a * |7(V’(s))I p )|7W’(>))I p = 0. (27) 

e_>0 J(Ofx (M«\Of) 

Equations (26) and (27) yield 

lim f (I a * |7(V’(s))| p )|7(V’(s))l P - [ (4 * |7(^i(s))| p )|7(^i(s))| p 
£_>U J R w J R w 

- f (i a *\j(Ms))\ p )h(Ms))\ p = o, 

J R N 

and hence, as £ —► 0, 

r e (^(s)) > r £ (^i(s)) + r e (^ 2 (s)) + o(i). 

By similar arguments, 


Teifois)) 

I I a (x-y)\'y(i/j 2 (s)(x)\ p \j(ip 2 {s)(y)\ p dxdy>o(l), 

Z P i(R»\O t )x(R w \O t ) 

and we finally conclude that 


r e{i>(s)) > r E (v>i(s)) + o(i) 


as £ —> 0. If we define 

4 >{(s)(x,y) 


ipi(s)(x,y) if x > 0 and y G O 26 
0 if x > 0 and y £ O 26 , 


we immediately see that r e (^>i(s)) > r E (-0{(s)), and ij}\ G <f> e because 0 < a < 
E Vo . Now [12, Proposition 3.4] implies that, as e —> 0, 


max r e (V>(s)) > E v +o(l), 

0 <s<T 


and this contradicts (24). 
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For a fixed e sufficiently small and for R 1, we consider a sequence 
{v R } n C Xe' 2 n rf' n HS(B{ 0, i?/e)) such that r^(u^) ->■ 0 is R/e)) 

as n —> +oo. The boundedness of {v R } n in Hq(B{ 0, R/e)) and the Sobolev 
embedding theorem imply that v R —► v R strongly in L q (B(0,R/e)) for any 
q < 2N/(N — 1). Since {v R } n is a Palais-Smale sequence, a standard argument 
shows that v R —► v R strongly in Hq(B( 0, R/e)). Hence the limit v R is a weak 
solution to the problem 


-Av R +m 2 v R = 0 in B (0, R ) 


with 


V e (y)v R (0, y) + (J Q * |u*(0,-)| p ) \v R (0,y)r 2 v R (0,y) + 

2p—l 

+ (2 p+l)([ Xel{v R f dy - l) XeV R {0,y ) 

VRJ'r /+ 

for y £ I w with |y| = R/e. 

Since v R £ X/ D Tj/ e n Hq(B( 0, R/e)), we deduce that both and 

{r e (u fi )}fl are uniformly bounded for e > 0 sufficiently small. Hence also 
{Qs(v r )}r is uniformly bounded for e > 0 sufficiently small. Now a Moser 
iteration scheme like [10, Theorem 3.2] yields that is bounded in L°° 

uniformly for s > 0 sufficiently small. Taking into account that is 

uniformly bounded in L°° and 

{I a * |u«(0,-)| p ) \v R (0 ,y)r 1 < hiy e +m)\v R (0,y)\ 

when |y| > 2 R/e, we can perform a comparison argument as in [10, Theorem 
5.1] and derive 

\v R {x,y)\ < Ce- m ^ x2+ \ y \ 2 ~ 2R o), 

We assume, without loss of generality, that {v r }r weakly converges to some v e 
in H as R —> +oo that solves 


9 »" 1 

-s£ (0 ' y) = - 


— Av e + m 2 v e = 0 in R^ +1 


(28) 


with 

f) 01 

--^(°,y) = ~V e (y)v £ {0,y) + (l a *\v e (0,-)\ p )\v e {0,y)\ p - 2 v e (0,y) + 

2p—l 

+ (2 p + l)([ Xelive) 2 dy - l) XeVe(0,y) (29) 
\J R« / + 

for y £R n . 
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6. Proof of the Theorem 1.1 


We can now collect all the results of the previous section to prove our main 
existence theorem. To begin with, Proposition 5.2 gives us a number £o > 0 
such that, for 0 < e < £q, the penalized functional r E possesses a critical point 
v e £ X£r\rf e . As in the proof of Proposition 3.3, we have v e £ H g >2 L q (R^ +1 ), 
and {i! e } is bounded L°°([ 0, +oo) x I*). By the results of Proposition 4.2, 

lim [ (|Vu E | 2 + V e \v e \ 2 ) dxdy = 0. 

£ ^°./R^ +1 \([0,+oo)x(.^ J ) £ ) 

It now follows that 

lim sup K(z,2 /)|=0 , 

(x,y)eB" +1 \([0,+oo)x(^ 2 3) c ) 

and as in the last step of the previous section we deduce an exponential decay 
of the trace u e away from \ (.# 2 ^) E : 

\u e (y)\ < Ci exp (-C 2 dist (y, („# 2/3 ) E )) . 

Taking e smaller, this estimate implies that Q s (v e ) = 0, and (28)-(29) are the 
local Neumann problem in the half space M. N corresponding to the nonlocal 
problem (4). The conclusion now follows by reversing the local realization of 
the operator y/— A + m 2 . Recalling (23) and all the scalings, we immediately 
deduce (5). This completes the proof. 
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